rd_mig7 <- RDestimate(mig ~ date | cntry, data = ESS10, cutpoint = 0, bw = 7, se.type = "HC1")
summary(rd_mig7)
rd_eco7 <- RDestimate(eco ~ date | cntry, data = ESS10, cutpoint = 0, bw = 7, se.type = "HC1")
summary(rd_eco7)
rd_cul7 <- RDestimate(cul ~ date | cntry, data = ESS10, cutpoint = 0, bw = 7, se.type = "HC1")
summary(rd_cul7)

rd_mig10 <- RDestimate(mig ~ date | cntry, data = ESS10, cutpoint = 0, bw = 10, se.type = "HC1")
summary(rd_mig10)
rd_eco10 <- RDestimate(eco ~ date | cntry, data = ESS10, cutpoint = 0, bw = 10, se.type = "HC1")
summary(rd_eco10)
rd_cul10 <- RDestimate(cul ~ date | cntry, data = ESS10, cutpoint = 0, bw = 10, se.type = "HC1")
summary(rd_cul10)

rd_mig14 <- RDestimate(mig ~ date | cntry, data = ESS10, cutpoint = 0, bw = 14, se.type = "HC1")
summary(rd_mig14)
rd_eco14 <- RDestimate(eco ~ date | cntry, data = ESS10, cutpoint = 0, bw = 14, se.type = "HC1")
summary(rd_eco14)
rd_cul14 <- RDestimate(cul ~ date | cntry, data = ESS10, cutpoint = 0, bw = 14, se.type = "HC1")
summary(rd_cul14)

rd_mig21 <- RDestimate(mig ~ date | cntry, data = ESS10, cutpoint = 0, bw = 21, se.type = "HC1")
summary(rd_mig21)
rd_eco21 <- RDestimate(eco ~ date | cntry, data = ESS10, cutpoint = 0, bw = 21, se.type = "HC1")
summary(rd_eco21)
rd_cul21 <- RDestimate(cul ~ date | cntry, data = ESS10, cutpoint = 0, bw = 21, se.type = "HC1")
summary(rd_cul21)

rd_mig28 <- RDestimate(mig ~ date | cntry, data = ESS10, cutpoint = 0, bw = 28, se.type = "HC1")
summary(rd_mig28)
rd_eco28 <- RDestimate(eco ~ date | cntry, data = ESS10, cutpoint = 0, bw = 28, se.type = "HC1")
summary(rd_eco28)
rd_cul28 <- RDestimate(cul ~ date | cntry, data = ESS10, cutpoint = 0, bw = 28, se.type = "HC1")
summary(rd_cul28)
